Fluctuations and correlations in a frustrated S = 1/2 square lattice with competing 
ferromagnetic and antiferromagnetic interactions: a /iSR study 
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Zero and longitudinal field /iSR measurements in Pb2VO(P04)2 and BaCdVO(PC>4)2 , two pro- 
totypes of the frustrated S — 1/2 square lattice model with competing ferromagnetic and antifer- 
romagnetic interactions, are presented. Both systems are observed to undergo a phase transition 
to a long-range magnetic order at TV — 3.46 K, for Pb2VO(P04)2 , and at TV — 0.99 K, for 
BaCdVO(P04)2 • In Pb2VO(PC>4)2 both the temperature dependence of the order parameter and 
the longitudinal relaxation rate above TV are consistent with a two-dimensional XY model. On the 
other hand, for BaCdVO(P04)2 , which lies very close to the magnetically disordered region of the 
phase diagram where a bond-nematic order was predicted, a peculiar logarithmic increase of the 
relaxation is observed above TV. In both systems a rather broad distribution of internal fields at 
the muon sites is noticed below TV- The origin of this distribution is discussed in the light of the 
/J.SK experiments already performed on S = 1/2 frustrated antiferromagnets on a square lattice. 

PACS numbers: 76.75.-K, 75.10.Jm, 75.40.Gb 
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I. INTRODUCTION 

Strongly correlated electron systems with competing 
interactions are known to show rather rich phase dia- 
grams, with crossovers or phase transitions which de- 
pend on the relative magnitude of the competing en- 
ergy scales^. In certain cases rather complex scenar- 
ios are observed together with the insurgence of novel 
phases of matter, a situation typically found in frus- 
trated magnets where, for instance, spin ice or spin liq- 
uid ground-states arise. 2 In these insulating systems ei- 
ther the geometry of the underlying spin lattice or the 
geometry of the interactions causes the suppression of 
the long-range magnetic order and the onset of exotic 
phases. For instance, for a frustrated S = 1/2 square lat- 
tice with competing ferromagnetic and antiferromagnetic 
interactions 3 (hereafter QFFMSL for short) it has been 
proposed that for exchange couplings yielding to the dis- 
appearance of long-range magnetic order, a nematic order 
could be present.— QFFMSL are characterized by a fer- 
romagnetic nearest neighbour (n.n.) exchange coupling 
J\, along the sides of the square spin lattice, compet- 
ing with the antiferromagnetic next nearest neighbour 
(n.n.n.) coupling J2, along the diagonals. For a frus- 
tration ratio —0.4 > J%/ J\ > —0.7 long-range magnetic 
order should disappear and a nematic order of spin bonds 
should arise^ In this phase, the individual spins do not 
carry any orientational order but the traceless tensor 
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with a, 7 = x,y,z and i,j running over spin lattice 
sites, does order and gives rise to stripy like correla- 
tions in the plane. This interesting theoretical result 
lacks of any experimental confirmation, basically because 
there is no material characterized by a frustration ratio 
J2/T1 ~ —0.5. Nevertheless, in the last years there has 



been a considerable effort to synthesize novel vanadates, 
containing V 4+ 5=1/2 ions, which have exchange cou- 
plings in this region of the J\ — J2 phase diagram. Several 
prototypes QFFMSL have been synthesized 5 , with ratios 
J2/J1 ranging from -1.8 for Pb 2 VO(P0 4 ) 2 to -0.9 for 
BaCdV0(P0 4 )2— , very close to the boundary between 
the long-range magnetic order and the nematic order. On 
the other hand, it has been recently pointed out by Tsir- 
lin and Rosner— that in these compounds the n.n. and 
n.n.n. exchange couplings along different sides and diag- 
onals of the square lattice, respectively, might be inequiv- 
alent due to the buckling and stretching of the magnetic 
layers. The difference between the exchange couplings 
appears to be less pronounced for BaCdVO(P04)2 than 
for Pb 2 VO(P0 4 ) 2 . 
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FIG. 1: Time evolution of the ZF/iSR asymmetry in 
Pb2VO(P04)2 at two representative temperatures, above and 
below TV. The solid lines are best fits according to Eq. 2 
(data at T = 3.7 K) and 3 (data at T = 3.37 K) in the text. 

In order to study the evolution of the local microscopic 
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properties of QFFMSL upon decreasing the ratio J2 / J\ 
towards the critical value J2/J1 — —0.7 and how they 
compare to the ones of frustrated 5 = 1/2 square lattice 
systems with competing antiferromagnetic interactions^' 9 
(QFAFSL, i.e. for J2/J1 > 0), we have performed 
zero field (ZF) and longitudinal field (LF) /zSR exper- 
iments on BaCdVO(P0 4 ) 2 and Pb 2 VO(P0 4 ) 2 powders. 
It was found that although both systems are character- 
ized by a magnetically ordered ground-state below Tn, 
significant differences with respect to the behaviour ob- 
served for QFAFSL are present. The correlated spin dy- 
namics in Pb2VO(PC>4)2 for T > Tm and the temper- 
ature dependence of the magnetic order parameter are 
both consistent with a two-dimensional (2D) XY model. 
On the other hand, in BaCdVO(P04)2 a rather pecu- 
liar behaviour is observed, with a logarithmic increase of 
the longitudinal muon relaxation rate (A) on cooling for 
T > Tn- This behaviour is reminiscent of the one ob- 
served in onc-dimensional systems and might suggest the 
onset of bond-nematic correlations above T/y. 
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FIG. 2: Power of the Fourier transform of the oscillating 
component of the ZF^iSR asymmetry for Pb 2 VO(P0 4 ) 2 , 
BaCdVO(P0 4 ) 2 and Li 2 VOSi0 4 , for T < T N . It is evident 
that while in Li 2 VOSi0 4 l — a well defined internal field is 
present, in Pb 2 VO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 a broad dis- 
tribution of internal fields is probed by the muons. 



II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

A polycrystalline sample of Pb2VO(P0 4 )2 was syn- 
thesized by solid state reaction technique using PbO 
(99.99%), (NH 4 )H 2 P0 4 (99.9%), and V0 2 (99.99%) as 
starting materials. In the first step, the intermediate 
compound Pb 2 P20y was prepared firing the stoichiomet- 
ric mixtures of PbO and (NH 4 )H 2 P0 4 in air at 750 °C 
for 24 hours. In the second step, the intermediate prod- 
uct was mixed with VO2 in appropriate molar ratio and 
heated for 48 hours at 680 °C in dynamic vacuum with 
one intermediate grinding and pelletization. Synthesis 
of BaCdVO(P0 4 )2 was done following the same proce- 
dure reported in Ref. [6j. Single phase materials were 
confirmed by x-ray diffraction performed with a STOE 
powder diffractometer (CuK Q radiation). Both materi- 
als, Pb 2 VO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 , contained a mi- 
nor (2-3%) fraction of unreacted diamagnetic phosphates 
Pb2P207 and BaCdP20y, respectively. These impurities 
are non-magnetic and, therefore, are irrelevant for the 
discussion to be presented in the following. 

ZF and LF/xSR measurements were performed on EMU 
and MUSR beam lines at ISIS pulsed muon facility, us- 
ing 29 MeV/c spin-polarized muons. The powders of 
Pb 2 VO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 were pressed and at- 
tached to a silver sample-holder, whose background con- 
tribution (B) to the muon asymmetry was determined 
from the slowly decaying part of the longitudinal polar- 
ization. We checked that there were no history dependent 
effects associated with a poor sample thermalization onto 
the dilution fridge cold plate. 

In Pb2VO(P0 4 ) 2 at temperatures above Tn = 3.46 ± 
0.01 K the decay of the muon polarization P M (t), either 
in zero (Fig. I) or in a longitudinal field, could be fitted 



by the expression 

f%(t) - Aexp(-(Xtf) + B . (2) 

In ZF, for T ^S> Tm, (3 was found to increase slightly 
above unity which indicates that at high temperature a 
non-negligible contribution to the relaxation arises from 
the field distribution generated by the dipolar interaction 
with the nuclei, which typically gives a Gaussian decay. 10 
On the other hand, when a LF of 1 kGauss was applied 
in order to quench the nuclear dipole contribution, (3 was 
found to decrease from unity to about 0.7 on approaching 
T N . 

Below Tjv ZF experiments showed the insurgence of os- 
cillations in the muon asymmetry (Fig. I) which clearly 
indicate the presence of a long-range magnetic order, 
which generates an internal field at the muon site B^. 
Accordingly one observes a precessional signal at a fre- 
quency uj^ = jfiBfj,, with 7^ the muon gyromagnetic ra- 
tio. Nevertheless, the initial decay of the muon polariza- 
tion can hardly be fitted by assuming a well defined field 
at the muon. On the other hand, one can fit P^it) over a 
broad temperature range below Tn with the expression 

P^t) = A^xpi-a^Joi^B^t) + A 2 exp(-(\tf) + B , 

where Jo{x) is the zero-th order Bessel function of the 
first kind, which characterizes the decay of the muon po- 
larization in the presence of a distribution of internal 
fields^ which is further broadened by the decay constant 
a. In fact, if one performs the Fourier transform of the 
oscillating component of P M (t) one finds quite a broad 
spectrum (Fig. 2). Upon decreasing the temperature 
below Tn one observes a loss of the initial asymmetry, 
possibly due to a further broadening of the field distri- 
bution. The temperature dependence of derived by 
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FIG. 3: Temperature dependence of the internal field at the 
muon in Pb2VO(P04)2 as derived from ZFpiSR experiments. 
In the inset the same data are reported as a function of the 
reduced temperature for T — * Tn . The solid line is the critical 
behaviour expected for a 2D XY model. 



fitting the ZF decay of the muon asymmetry with Eq. 3 
is reported in Fig. 3. In Fig. 4 the temperature depen- 
dence of the longitudinal relaxation rate A is reported. A 
clear divergence of the relaxation rate is observed at Tjv. 

Also for BaCdVO(PC>4)2 at temperatures above Tjv = 
0.99 ± 0.01 K the decay of the muon polarization can be 
nicely fitted with Eq. 2, with an exponent decreasing 
from j3 = 1 for T ~ 10Kto/3 = 0.6on approaching 
Tjv (Fig. 5). On the other hand, although in ZF be- 
low Tjv oscillations in P M (£) are observed, confirming the 
onset of a long-range magnetic order—, the decay of the 
muon polarization can hardly be fitted by Eq. 3 and a 
much broader distribution of internal fields is present. In 
fact, by looking at the Fourier transform of the oscillating 
component of the muon asymmetry (Fig. 2) one clearly 
notices that a significant distribution of internal fields is 
present, much broader than in QFAFSLi^. Still, it is pos- 
sible to fit the slowly decaying non-oscillating component 
of (t) with a stretched exponential in order to derive 
the temperature dependence of A over all the explored 
temperature range. Again a peak is present at Tjv (Fig. 
6) . However it is noticed that for T > Tjv the increase of 
A on cooling is rather different from the one observed in 
Pb 2 VO(P0 4 ) 2 ■ 



III. DISCUSSION 

First we shall address the temperature dependence of 
the muon longitudinal relaxation rate and the tempera- 
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FIG. 4: Temperature dependence of the muon longitudinal 
relaxation rate A in Ph)2VO(P04)2 . Above Tjv A was derived 
in a field of 1000 Gauss, while below Tjv it was estimated from 
the non-oscillating component of the muon asymmetry. The 
solid line shows the behaviour expected for a 2D XY model. 
In the inset the same data are reported for T — > Tat as a 
function of \J J C /(T — Tn), together with the data derived for 
Sr2Cu02Cl2, for a LF H — 800 Gauss along the c axes. J c = 
10.4 K for Pb 2 VO(P0 4 )2 and J c = 1450 K for Sr 2 Cu0 2 Cl2.— 
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FIG. 5: Time evolution of the ZF/iSR asymmetry in 
BaCdVO(P04)2 at three representative temperatures above 
Tjv. The solid lines are best fits according to Eq. 2. 



ture dependence of i? M in Pb2VO(PC>4)2 ■ The temper- 
ature dependence of corresponds to the one of the 
order parameter. In fact, in zero-field the local field at 
the muon can be written as = A < S >i£, where A 
is the hyperfine coupling tensor describing the coupling 
of the muon with the surrounding V 4+ S = 1/2 spins 
and < S > the corresponding spontaneous spin polar- 
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ization. In spite of the distribution of local fields evi- 
denced in Fig. 2, which yields some uncertainty in the 
estimate of B^, the basic trend of B^iT) is reminiscent 
of the one observed for Li 2 VOSi0 4 QFAFSLJ^ In this 
latter compound it has been observed that for T — > Tjv 
the order parameter increases according to the critical 
power law B^ oc (1 — T/Tv)' 30 , with a critical exponent 
j3 c ~ 0.235, the one predicted for a 2D XY system*^ 
If one reports B^ data for Pb2VO(P04)2in the prox- 
imity of Tjv as a function of the reduced temperature 
(1 — T/Tv) (see the inset of Fig. 3), one realizes that 
also for Pb 2 VO(P0 4 ) 2 & ^ 0.235. 
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FIG. 6: Temperature dependence of the zero-field muon lon- 
gitudinal relaxation rate A in BaCdVO(P04)2 . In the inset 
the same data are reported vs. 1/T in a linear- log scale in 
order to evidence the logarithmic increase of A above Tjv. 

The temperature dependence of the muon longitudinal 
relaxation rate above Tjv originates from the one of the 
in-plane correlation length. In fact, if we neglect the nu- 
clear contribution to the relaxation which is temperature 
independent and in anyway, quenched by the application 
of a longitudinal field, the critical increase of A on ap- 
proaching Tjv is of dynamical nature and associated with 
spin-lattice relaxation processes^ Then one can write 

A =^ElAf^ Q (?,^) , (4) 

q.a 

where S aa (q,LO^) are the components of the dynamical 
structure factor at the muon Larmor frequency and 
is the form factor describing the hyperfine coupling of 
the spin excitations at wave-vector q with the nuclei. 
By resorting to scaling arguments, which are expected 
to be valid for T — > Tv, one finds that for a 2D spin 
system A cx £ z , with £ the in-plane correlation length 



and z the dynamical scaling exponent. 14 Considering that 
Pb 2 VO(P04) 2 lyes in the sector of the J\ — J 2 phase di- 
agram where the antiferromagnetic J 2 coupling is larger, 
one could at first take for £ the temperature dependence 
expected for a S = 1/2 2D Hcisenberg antiferromagnet 
on a square lattice. Namely, £ oc exp(2irp s /T)l^, with 
p s the spin-stiffness constant, which is possibly reduced 
with respect to the value p s = 1.15J/2tt expected for 
a non-frustrated systems. Accordingly one should ob- 
serve that A oc exp(2zirp s /T). However, we have found 
that this expression can fit the data only over a lim- 
ited temperature range, no matter which value one takes 
for p s . On the other hand, the temperature dependence 
of the order parameter below Tv suggests considering 
for £ the form predicted for a 2D XY modet^, namely 
£ oc exp[B X y I '(T — Tv)] 1 / 2 , where B X y is of the order 
of magnitude of the exchange coupling constant. Then, 
one should find 

A = CX e tfB XY /{T-T N )]V* (5) 

Indeed, Pb 2 VO(P04) 2 data can be nicely fitted with 
Eq. 5 over almost all the explored temperature range 
(Fig. 4). It is worth comparing the behaviour ob- 
served for Pb 2 VO(PC>4) 2 with the one found for a pro- 
totype of the S = 1/2 2D XY model, Sr 2 Cu0 2 Cl 2 Ji In 
Fig. 4 we compare the temperature dependence of A for 
Pb 2 VO(PC>4) 2 and for Sr 2 Cu0 2 Cl 2 , reported as a func- 
tion of J C /{T — Tjv), with J c = y/ J\ + jf an effective 
exchange coupling. One notices that for both systems 
Eq. 5 is followed very well. 

Now we shall address the temperature de- 
pendence of the order parameter describing the 
BaCdVO(P04) 2 ground-state. A broad distribution of 
local fields at the muon site, even larger than the one 
found for Pb 2 VO(P0 4 ) 2 , is noticed (Fig. 4). This 
field distribution should not be associated with sample 
inhomogeneities since the phase transitions detected by 
A or specific heat measurements 6 appear rather sharp. 
On the other hand, it should be mentioned that in 
QFAFSL as Li 2 VOSi04 , where the buckling along the 
a-axis is not present and the unit cell is tetragonal, such 
a broad distribution is absen t 12 ! 18 (Fig. 4). Thus, at 
first sight one could think that the different structure 
could yield a distribution of muon sites, leading to a 
distribution of hyperfine coupling constants which could 
justify the broadening of the Fourier transform spectra. 
However, it is rather difficult to associate such a broad 
distribution with the presence of a discrete number 
of p + sites and moreover the distribution is larger in 
BaCdVO(P04) 2 where the ab plane is less distorted. 7 
Furthermore, even if a priori the presence of several 
incquivalcnt muon sites cannot be excluded, we mention 
that in oxides with superstructures and much larger 
unit cells, as Bi 2 Sr 2 YCu 2 0g-^ and Sri4Cu 2 404i— for 
instance, just up to two muon sites are observed. 

If the significant broadening in the Fourier transform is 
not due to a distribution of inequivalent probes it means 
that < S > varies from lattice cell to lattice cell, namely 
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either the underlying magnetic lattice is not commen- 
surate with the crystal lattice structure or domains are 
present at the microscopic level. In this respect it is 
instructive to compare our results to neutron scatter- 
ing onesi^l Although neutron scattering experiments in 
Pb2VO(PC>4)2 indicate the presence of a columnar order, 
a broad background distribution of the energy integrated 
dynamical structure factor is noticed in q-space, even well 
below Tjvi^i The origin of this broad distribution could 
suggest the presence of small size domain like structures. 
Indeed Pb 2 VO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 are in a sec- 
tor of the J\ — J 2 phase diagram where a two-fold degen- 
erate magnetic ground-state is present^ and one could 
speculate that under certain conditions the coexistence of 
domains of both phases, characterized by different mag- 
netic wave vectors, could be observed. Although, we can- 
not draw a clear conclusion on the origin of the broaden- 
ing of the spectra, it is clear that in BaCdVO(P04) 2 and 
Pb 2 VO(P04) 2 a larger degree of structural and possibly 
of magnetic disorder is present with respect to QFAFSL 
as Li 2 VOSi0 4 • 

In BaCdVO(P04) 2 the application of a moderate mag- 
netic field is observed to affect significantly the spectrum 
of the excitations and the specific heati£ Therefore, in 
order to probe the intrinsic low-energy excitations of this 
system we have performed zero-field relaxation experi- 
ments. Under these conditions an additional contribu- 
tion to the decay of the muon polarization, associated 
with the dipolar interaction with the nuclei, might be 
present. Nevertheless, the functional form of -P M (i) sug- 
gests that in BaCdVO(P0 4 ) 2 this contribution is less sig- 
nificant than in Pb 2 VO(P0 4 ) 2 , where the nuclear con- 
tribution to the relaxation is already minor for T — > Tn- 
In fact, in this latter compound A is observed to diverge 
for T — > Tjv according to Eq.5 also in ZF. 

On the other hand, the temperature dependence of 
the zero-field relaxation rate in BaCdVO(PC>4) 2 shows 
a rather peculiar behaviour, not observed in 
Pb 2 VO(P0 4 ) 2 . At high temperature, for T > 4 
K, A is temperature independent, as expected for a 
non-correlated spin system. On decreasing the temper- 
ature towards Tjv, one would expect for a 2D system, 
on the basis of Eq. 4, an exponential growth of A. 
Remarkably A is found to increase logarithmically, i.e. 
A oc ln(X/T) (Fig. 6), a trend which is rather charac- 
teristic of one-dimensional systems^ This behaviour 



should not be related to the filtering of the critical 
fluctuations by the form factor in Eq. 4^ In fact, if this 
was the case one should expect the cancellation of the 
magnetic field at the muon site below Tjv, at variance 
with the experimental findings (Fig. 2). On the other 
hand, since BaCdVO(P0 4 ) 2 should lie very close to the 
boundary with the bond-nematic phased, one could 
expect that nematic correlations of stripy character 
arise above Tn and yield a logarithmic increase of A 
with decreasing temperature. This occurs until the 
XY anisotropy and/or the interlayer coupling cause the 
crossover to a three-dimensional long-range order at a 
finite temperature resulting in an abrupt increase of A 
very close to Tn (see Fig. 6). 



IV. CONCLUSIONS 

In this manuscript we have presented new ^SR 
measurements in two prototypes of QFFMSL: 
Pb 2 VO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 . In both com- 
pounds a broad distribution of local fields at the muon 
site is evidenced and is tentatively associated either 
with an incommensurate magnetic order or with the 
formation of mesoscopic domains. In Pb 2 VO(P04) 2 the 
overall temperature dependence of the longitudinal 
relaxation rate and of the order parameter below Tjv 
are both consistent with a 2D XY model. On the other 
hand, in BaCdVO(P04) 2 a rather peculiar behaviour of 
the longitudinal relaxation rate was evidenced, which 
can be hardly associated with conventional 2D spin 
correlations. The logarithmic increase of A on cooling 
suggests the onset of one-dimensional correlations which 
might be related to the insurgence of novel type of 
correlations in this QFFMSL. 
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